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ABSTRACT 

Crystal oscillators are now available with frequency stabilities 
of a few parts in i0°° or ion. One application of such stability 
is in navigation systems, where geographic position may be determined 
by phase-comparing a received frequency-stable signal with a local 
frequency-sStable signal. Errors may be introduced, however, by the 
effects of static position and of motion on the frequency stability 
of the oscillators employed. A testing program including static posi- 
tion, swinging motion, vibration, and shock tests was devised to 


determine the extent of such effects. The program is described and 


the results are presented in both tabular and graphic form. 
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i, Introduction. 

Due to improvements in crystal-unit and oscillator design and 
construction during the past ten years, it is now not unusual to 
speak of crystal oscillators with the capability of maintaining their 


11 
or 10 , Oscil- 


frequency outputs constant to a few parts in 10! 
lators of this sort have been employed in very low frequency naviga- 
tion systems, where the phenomenon of long-range, phase-stable VLF 
transmissions coupled with frequency-stable signals at the trans- 
mitters and aboard ship, permits a geographic position to be determined 
by a process of phase-comparison. 

The accuracy of such a navigation system depends upon the 
frequency stability of the oscillators, It is recognized that en- 
vironmental conditions such as temperature, vibration, and shock have 
a direct effect on the operation of crystal oscillators and, therefore, 
every effort is generally made to provide the operating conditions 
as specified by the manufacturer of the equipment. 

In addition to the problems of temperature, vibration, and shock 
it was suSpected that motion such as that created by the rolling 
and pitching of a ship or buoy might introduce instability. 

It was the purpose of this investigation to determine the effects 
of such motion on the frequency stability of certain crystal oscil- 


lators as well as the effects of static position, vibration, and shock, 





2. Oscillators Tested, 
Tests were conducted with the following oscillators: 
RF Oscillator 0-76A/U 
Western Electric D-175730-L2 Frequency Standard 
Serial No. 43 


Sulzer Laboratories Model D5 Oscillator 
Serial No. 25 


Sulzer Laboratories Model D5 Oscillator 
Serial No, 19 


Sulzer Laboratories Model D5 Oscillator 
Serial No. 22 


Sulzer Laboratories Model 2.5 Frequency Standard 
Serial No. 256 


A brief description of each type of oscillator is given in the 
following paragraphs, Further information may be obtained from the 
equipment manuals furnished by the manufacturers. [1, 6, 7\. 

The Western Electric frequency standard produces a frequency 
of 100 kilocycles which, according to the manufacturer, can be 
maintained constant to within approximately one part in 10° per day. 
The crystal, a GI-cut quartz plate with a very low temperature co- 
efficient, is maintained at a constant temperature by a two-stage 
oven and a bridge oscillator circuit. The oscillator weighs 90 
pounds; its dimensions are 19" x 20" x 10", 

The Sulzer Laboratories Model D5 oscillator is a solid-state 
unit with a five-megacycle output and a rated stability of five 
parts in Ae per day. The crystal is an AT-cut quartz plate, vibrating 
in a fifth-overtone, thickness-shear mode, Temperature level is 


maintained by a two-stage oven. The oscillator weighs approximately 


are 





three and one-half pounds; its dimensions are 4%" x 4%" x 6 7/8". 
The third type of oscillator tested, a Sulzer Laboratories 

Model 2.5 frequency standard, is also a solid-state unit. It has 
frequency outputs of two and one-half megacycles, one megacycle, and 
100 kilocycles, Stability as stated by the manufacturer is one part 
in 10°” or better per day. The crystal-unit is a two and one-half 
megacycle military-type designed by Warner of Bell Telephone Labora- 
tories, It is mounted in a proportionally-controlled double oven. 
The oscillator weighs approximately six and one-half pounds; its 


dimensions are 4%'' x 44" x 11 3/4", 





oe Test Program, 

Purpose and Procedure 

The purpose of the testing program was to reveal the effects 
of static position, swinging motion, vibration, and shock on the 
frequency stability of the crystal oscillators described above. To 
accomplish this goal, each of the oscillators was subjected to certain 
conditions and combinations of position and motion while a constant 
observation was made of its frequency output. 

The frequency standard employed as a reference for the tests 
was a Varian Model V-4700 Rubidium Vapor Frequency Standard with a 
long term frequency stability of 5 x tas any one-year period 
(standard deviation). To insure the reliability of the standard 
during the test period, its output was continually phase-compared 
against the 20 kilocycle signal of WWVL, This comparison was made 
with a Textran Series 599 VLF Tracking Receiver. The stability 
during the period from 23 July 1965 to 1 November 1965 was 2.3 x Gees 
(standard deviation). 

The output of the test oscillator was phase-compared against 
the Rubidium standard and the output of the phase comparator was 
recorded on a strip chart recorder, The basic test system is shown 
in Figure 1. Examination and processing of the recorded data 


produced the results presented in Section 5, 





Circuits and Equi pment 


While the test circuit was basically that shown in Figure l, 
the laboratory hook-up was more complicated due to the simultaneous 
testing of as many as four oscillators. All phase comparisons were 
made at 100 kilocycles, necessitating the use of frequency dividers 
in the case of the Sulzer D5's with their five-megacycle outputs. 

The phase comparators, frequency dividers, and oscillators also 
required separate power supplies. Figure 2 shows an extension of 
the basic test system for the case of four oscillators under test. 

A thermistor-bridge circuit and a chart recorder were employed 
to maintain a constant record of room temperature for the duration of 
the tests, 

Swinging motion was provided by a specially constructed swinging 
platform shown in Figures 3 and 4, The swing was driven by a three- 
quarter horsepower electric motor coupled to a hydraulic transmission. 
. The transmission permitted adjustment of swing frequency and the 
slotted drive wheel permitted adjustment of swing amplitude, 

Vibration was accomplished with a Goodmans Model 390 A Vibration 
Generator. A low-frequency function generator and a power amplifier 
were used to supply the necessary input. The shake table with the 
three Sulzer D5 oscillators mounted on it is shown in Figure 5, 

Figure 6 illustrates the equipment used to provide a periodic 
lift-and-release motion for the shock testing described in a later 
paragraph, The equipment consisted of a direct-current motor with 


a cam attached to its shaft, and geared for high-torque, slow-speed 
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BASIC TEST. SYSTEM EXTENDED FOR 
SIMULTANEOUS TESTING OF FOUR SULZER OSCILLATORS 
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FIGURE 3 
SWING PLATFORM, VIEW 1 
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SWING PLATFORM, VIEW 2 











FIGURE 5 
SULZER D5 OSCILLATORS MOUNTED ON SHAKE TABLE 
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operation. The cam periodically lifted and released a hinged bar 
from which was suspended the mounting platform for the oscillators. 

Description of Tests 

Static Tests, The static tests consisted of operating the test 
oscillator for at least 48 hours on each of its six sides. The 
oscillator was operated first in its normal position, then rotated. 
counterclockwise through the side, upside down, and side positions. 
It was then operated on each of the remaining two sides, After 48 
hours in each of these positions, the oscillator was returned to its 
normal position for an additional 48 hours of testing. These tests 
were conducted on all of the oscillators. Table I lists the various 


Static tests and their descriptions, 


TABLE I 


STATIC TEST DESCRIPTIONS 


re or a TT eT 


Test Oscillator Position 





Static One Oscillator in normal operating position 

Static Two Oscillator tipped onto its side 

Static Three Oscillator tipped upside down 

Static Four Oscillator tipped onto its side 

Static Five Oscillator tipped forward onto its face 

Static Six Oscillator tipped backward onto its back 








Swing Tests. These tests were also conducted over a basic 
48 hour test period. The test oscillator or oscillators were 
fastened to the swing platform and the swing was then placed in 
motion, The Western Electric 0-76A/U was swung at various ampli- 
tudes and frequencies in one position only, the normal position, 
with the front panel of the oscillator parallel to the leading edge 
of the swing platform and perpendicular to the line of motion, The 
four Sulzer oscillators were swung simultaneously, and in two posi- 
tions, the first being identical to that described for the 0-76A/U. 
For the second position, the oscillators were rotated 90° in the 
plane of the platform, The period and amplitude of the various tests 


are given in Tables II and III. 


TABLE II 


SWING TEST DESCRIPTIONS, 
WESTERN ELECTRIC OSCILLATOR 


ee eS PP A aS AE 





Test Amplitude* Period 
fe) fe) 
Swing One 10°; 10° 30 seconds 
Swing Two 10,, 105 10 seconds 
Swing Three 20,, 29 30 seconds 
Swing Four Jo: 5 30 seconds 
Swing Five 259 b, 7 10 seconds 
Swing Six FO gg, tO 20 seconds 


a s 


*10°, ito represents a symmetric swing of 10° off the vertical 
in each direction. 


10 





TABLE III 


SWING TEST DESCRIPTIONS, 
SULZER OSCILLATORS 


—— 


Test Amplitude* Period 





First Position: 


Swing One Ss 10° 30 seconds 
Swing Two 5 aw 10 seconds 
Swing Three 24°, 30. 30 seconds 
Swing Four 24°, 30 10 seconds 


Second Position: 


Swing Five a 10° 30 seconds 
Swing Six 36: 10, 10 seconds 
Swing Seven 24 0s 305 30 seconds 
Swing Eight 24°, 30 10 seconds 


O O ; : O ; 
*5°, 10 represents an aSymmetric swing of 5. off vertical in 


one direction and 10° off vertical in the opposite direction 
Vibration Test, This test was conducted only with the Sulzer 
D5 oscillators. The osScillators were mounted on a one-quarter inch 
thick aluminum plate with a diameter of twelve inches which was then 
attached to the vibration generator (see Figure 5). They were vi- 
brated simultaneously at a frequency of 20 cycles per second with a 
force of approximately one g for a period of 12 hours. After a 
12 hour quiet period, the oscillators were once again subjected to 
the same vibration condition for a period of ten hours. 
Shock Test. The shock, or impact, test was also conducted only 
with the Sulzer D5 oscillators. The test consisted of a free fall 


from a height of three-tenths of an inch onto asphalt tile laid on 


yk 


concrete (see Figure 6). As in the vibration test, all three 
oscillators were tested simultaneously. Impacting occurred at a 
rate of once every ten seconds and was conducted continuously for 
12 hours, 

Environment. All tests were conducted in an open laboratory 
with room temperature normally ranging daily from 21°C to 20mer Over 
the 43 day test period the mean temperature was 22°C with a low of 
18°C and a high of 29°C. 

Data furnished by Sulzer Laboratories with respect to effects 
of temperature change on frequency stability of the Sulzer oscillators 
is shown in Table IV. Similar data for the Western Electric oscil- 
lator was not available, 

TABLE IV 


EFFECTS OF TEMPERATURE ON 
FREQUENCY STABILITY OF SULZER OSCILLATORS* 








Oscillator Temperature Maximum Af** 
Change F 
-10 
Model 2.5 ambient temperature + 2 x 10 
Frequency Standard 20°C changed + 15 Sc 
Model D5 ambient temperature + 1x nome 
-54°C vs + 25°C 
Model D5 ambient temper ature ro XO at 
-18°C vs + 25°C 
Model D5 ambient temperature + 5 x 10 ae 
+60°C vs + 25°C 
Model D5 ambient temperature i 10" 


+65°C vs + 25°C 
TS A 
*Information furnished by Sulzer Laboratories 
*kDefined on Page |4 


--— 


4, Data Acquisition and Processing. 

Acquisition 

During the, course of each test, the outputs of the various phase 
comparators were fed to strip chart recorders which provided a 
permanent data record for processing and analysis at a later date. 
Temperature data was also permanently recorded. The phase-comparison | 


data was in the form shown in Figure 7. 
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SAMPLE OF RECORDED PHASE COMPARATOR OUTPUT 
The output of a phase comparator fea f@ett the difference between 
“the frequencies of the two input signals. If the frequencies of the 
signals were equal, the phase difference of the signals was constant, 
and the output of the phase comparator was constant, If the fre- 
quencies were not equal, however, the signals passed in and out of 
phase with each other periodically and the cyclic change in phase 
difference produced a similar variation in the phase comparator output. 
By recording this Eaite output on a strip-chart moving at a regulated 
speed, the ite of phase shift, cycles/unit ae determined. 


For example, as shown in Figure 7, during the period from 2100 to 2200 


13 





the phase comparator passed through 1.05 cycles; the rate of phase 
shift during this period was 1.05 cycles/hour. 

All phase comparison was done at a frequency of 100 kilocycles. 
Therefore, during the one-hour period, the standard signal went 
through 3600 x 100,000 or 3.6 x 10° cycles while the test signal 
went through 3.6 x 10° - 1.05 cycles. 

The difference between the frequencies, expressed as a fraction, 


is called the fractional frequency error or the frequency offset. 


It is written as Af where 
F 
f = frequency of test signal 


F = frequency of standard signal 
Af=f-F 
The frequency offset for the above example was 


a ee rT See 


. 3.6 x 10 
It is customary to refer to the amount of phase shift in terms 
of time rather than frequency. At 100 kilocycles, the time for one 
cycle is ten micro-seconds and 1.05 cycles of phase shift per hour 
is equivalent to 10.5 micro-seconds of phase shift per hour. The 
frequency offset may then also be computed as 


Af = At = -10.5 x 107° = -2.9 x 10. 


F £ 3600 


9 


where T = period of observation in seconds 
At = phase shift in seconds occurring during T 


Processing 


The recorded data for each test was examined in the manner out- 


lined above. The number of micro-seconds of phase shift per hour 


14 





was noted for each hour of each test period. These figures were 
used to compute hourly values of frequency offset, Af/F, which 
were then plotted versus time. A sample plot is shown in Figure 8. 
The positive values of frequency offset indicate that the frequency 
of the test oscillator was above that of the standard, The upward 
drift of the points represents an increase in frequency of the test 
oscillator. This is a normal condition, caused by aging of the 
crystal. 

Next, a digital computer was employed to fit a straight line 
to the plotted points and to measure the dispersion of the points 
about this fitted line by computing their standard deviation. The 
program used for these computations is given in the Appendix. The 
Slope of the fitted line represented the frequency drift rate of the 
test oscillator and the standard deviation of the points was an 
indicator of the oscillator's frequency stability, with lower values 
of standard deviation indicating better stability. Figure 9 shows 
the sample plot presented in Figure 8 with a straight line fitted 
to the points. The choice of a one-hour interval between the points 
was purely arbitrary. Figures 10 and 11 present the same data plotted 
at intervals of one-half hour and two hours respectively. It may be 
noted that the smallest variation of Af/F plotted in Figure 9 
is / x TOR eae in Figure 10, 14 x or and in Figure 11, 3.5 x ome 
This is because the smallest value of Af/F is inversely propor- 
tional to the length of the observation interval. Using the 40-line 


chart paper shown in Figure 7, with the recorder calibrated for ten 
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micro-seconds of phase shift for full scale deflection, the smallest 
increment of phase shift readable is 0.25 micro-second. Therefore, 
for a one-hour observation interval, the smallest value of Af/F 
(other than zero) is 


6 et 


Af = At = 0,25 x 10. = 7x10 
F sh 3600 


For the purposes of this investigation, the choice of a one-hour 
observation interval proved to be quite satisfactory. Additional 
accuracy could have been obtained, however, by using recording equip- 
ment with 100-line paper, permitting phase shift to be read to 0.1 
micro-second, 

In summary, the data processing was accomplished in three steps. 
First, the strip charts were examined for the amount of phase shift 
per hour. Second, this data was converted to values of frequency 
offset and plotted against time. Third, to the plotted points was 
fitted a straight line, whose slope was thus obtained, and a measure 


of the variation of the points about the line was computed. 
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a. Test Results. 

The results of all tests are presented in Tables V through 
IX and Figures 12 through 37. The complete testing program of 
Static, Swinging motion, vibration, and shock tests was conducted 
with only the three Sulzer Model D5 oscillators. The two ee 
standards were not subjected to vibration and shock. 

Prior to the arrival of the Sulzer equipment, static and 
Swinging motion tests were begun with two Western Electric fre- 
quency standards. One of these was the oscillator mentioned in 
Section 2, a model 0-76A/U, and the second was an earlier model, 
076/U. It was soon determined that the drift rate and frequency 
stability of the older 0-76/U were affected by fluctuations in room 
temperature, The combined weight of the two oscillators (approxi- 
mately 180 pounds) also proved to be too great a load for the swing 
platform and prevented simultaneous swing testing. A circuit 
failure in the 0-76/U, coupled with the problems of its poor stabi- 
lity and load on the swing, caused testing with this particular 
obcillator to be discontinued. The 0-76A/U, however, was unaffected 
by changes in room temperature and operated properly during the test 
period. 

The results of the static and swing tests for this oscillator 
are given in Table V and Figures 12 through 15. In some instances 
the movement of the oscillator from one static position to another 
created a temporurzy frequency disturbance and several hours were 


required for the drift rate to return to its former state. 
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During this recovery period the change in frequency offset was 
relatively rapid and therefore the total test period was divided 

into sections to permit better fitting of straight lines to the 

data points. The three highest values of drift rate plotted in 
Figure 12, and the three highest values of standard deviation plotted 
in Figure 13, occurred during the recovery periods. For reference 
purposes the oscillator was operated in the Static One or normal 
position at the beginning and end of both the static test series 

and swing test series. The "motor only" entry in Table V refers to 

a period when the swing drive was disconnected and the swing platform 
was at rest, but the drive motor was running. This particular test 
was conducted to check for possible effects of vibration generated by 
the motor and hydraulic transmission. 

After a three-week warm-up period, the new Sulzer equipment was 
tested, It was recognized that this was an inadequate warm-up time, 
but the late arrival of the equipment and the pressure of time made 
it mandatory to begin the tests. The results of these tests are given 
in Tables VI through IX and Figures 16 through 37. The room tempera- 
ture variations had no noticeable effect on the drift rate or fre- 


quency stability of these oscillators. 
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TABLE V 


TEST RESULTS FOR 
WESTERN ELECTRIC 0-76A/U OSCILLATOR 
SERIAL NO, 43 


Test Condition Length Drift Rate 9 Standard 
(Hours) Af/F (pp 10 ) Deviation 
hour (x 10°) 
Static 1 normal 48 0.110 O32 
Static 2 (a) side 24 0.473 ees 
Seatic 2 (b) side 24 0.141 OFS 
Static 3 upside down 48 0.143 OTe 
Static 4 side 40 Onl 72 O55 
Static 5 (a) front 2 0.835 Vek2 
Static 5 (b) front 36 0.053 0.44 
Static 6 (a) back 6 1.494 Pot 
Static 6 (b) back 20 0.237 0.49 
Stacic, 6 (c) back 22 0.144 0.43 
Static l normal 48 0.102 0.82 
Static’ | normal 192 02037 Ole 7; 
Static 1 normal 158 0.025 0.41 
Swing 1 10°, 10°, 30 sec.* 48 0.133 0.80 
Swing 2 10°, 10°, 10 sec. 48 0.033 0.80 
Swing 3 20_, 20 , 30 sec. 48 0.031 On7Z 
Swing 4 5°, 5°, 30 sec. 48 -0,003 0.86 
Swing 5 35: 5 3 10 sec. 48 0,047 Lez 
Swing 6 10 , 10 , 20 sec. 48 -0.037 0.81 
Static 1 (motor only) normal 48 0.039 0.67 
Static | normal 48 O2127 OFZ 6 
Static 1 normal 48 0.033 0.60 
Static l normal 97 0.025 0.73 
SS ae 
nae O O 
10°, 10 , 30 sec. represents a symmetric swing of 10 off 


the vertical in each direction with a 30-second period. 
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TABLE VI 


TEST RESULTS FOR 
SULZER MODEL 2.5 FREQUENCY STANDARD 
SERIAL NO. 256 











Test Condition Length Drift Rate 11 Standard 
(Hours) Af/F (pp 10 _) Deviation 
hour (x 10°) 
Static 1 normal 65 2.96 5.04 
Static 2 side 49 7a EN. 6.16 
Static 3 upside down 49 2205 5.41 
Static 4 side 48 Oe, 4.1/7 
Static 5 front 48 t 60 Sle (a 
Static 6 back 48 Lis 1 
Static 1 normal 68 1.47 ee 18, 
Swing 1 5°, 10°, 30 sec.* 49 1.26 7.95 
Swing 2 36> 10°; 10 sec, 50 lee 9.40 
Swing 3 24 0» 30); 30 sec. 47 1,04 207 
Swing 4 24 0s 305: 10 sec, 47 1,09 6.97 
Swing 5 ns 10°; 30 sec. 48 8) 2)! 559) 
Swing 6 353 10.5 10 sec. 48 Or e)e Je) 
Swing 7 24 , 30 _, 30 sec. 48 OF FAG. 4.46 
Swing 8 24°, 30°, 10 sec. 48 OMT 4.83 
Static 1 normal 48 Oral Sh [pe 
O O ; O 
*5 , 10 , 30 sec. represents an asymmetric swing of 5 off 


vertical in one direction and 10 off vertical in the opposite 
direction with a period of 30 seconds. 
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TABLE VII 


TEST RESULTS FOR 
SULZER MODEL D5 OSCILLATOR 
| SERIAL NO. 25 





Test Condition Length Drift Rate il Standard 
(hours) Af/F (pp 10 -) Deviation 
hour (x 10°) 
Beatic 1 normal 65 -4,14 8.46 
Static 2 side 49 -3.61 (26 
Static 3 upside down 49 <3 4h3 Oral / 
Static 4 side 48 201 ie 18, 
Static 5 Front 48 -2.51 5.93 
Static 6 back 48 -2.26 O. ae 
Static l normal 68 =2,31 5.66 
Swing 1 5, 10), 30 sec.* 49 -2,73 15.06 
Swing 2 - 10., 10 sec 50 -2,25 12.63 
Swing 3 24 30); 30 sec, 47 -1.90 12.63 
Swing 4 240 30, 10 sec, 47 -1.47 4,28 
Swing 5 5 ., 10 _, 30 sec 48 -1,41 6.64 
Swing 6 5°, 10°, 10 sec 48 -1,42 9,28 
Swing 7 24 4s 300; 30 sec, 48 -1,61 ll. 2 
Swing 8 24 , 30 , 10 sec. 48 -0.68 9.30 
Static 1 normal 48 -0,81 1G322 
Vibrate 1 20 cps, 1 g iy =5;02 Se) 
Static 1 normal 12 -3,59 eZ 
Vibrate 2 20 cps, 1 g 10 -2,.42 9,64 
Static l normal 18 -3,01 8.45 
Static 1 normal ly -1.01 1209 
Shock 1 OFS" tall 12 +0,71 9 el 
Static 1 normal 60 -1.63 4,89 
45° 10°, 30 sec. represents an asymmetric swing of 5° off 


vertical in one direction and 10 off vertical in the opposite 
direction with a period of 30 seconds, 
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TABLE VIII 


TEST RESULTS FOR 
SULZER MODEL D5 OSCILLATOR 
SERIAL NO, 19 





Ce le 


Test Condition Length Pritt Rate, Standard 

(hours) Af/F (pp 10) Deviatign 

hour (x 10s) 
Static 1 normal 65 -6.29 8.46 
Static 2 side 49 -4,82 Seas 
Static 3 upside down 49 ~3.96 S250 
Static 4 side 48 -3.52 S.52 
Stacie 5 front 48 -3.02 8,86 
Static 6 back 48 -3,28 8.98 
Static l normal 68 -2.86 6.32 
Swing 1 5°, 10°, 30 sec.* 49 -3.31 8.75 
Swing 2 5°, 10°, 10 sec. 50 -2.68 7.68 
Swing 3 24 _, 30 _, 30 sec. 47 -2.46 4,74 
Swing 4 24°, 30°, 10 sec. 47 2a 4.12 
Swing 5 5, 10°, 30 sec. 48 -2.03 6.06 
Swing 6 5°, 10°, 10 sec. 48 SUN: 4.67 
Swing 7 24°, 30°, 30 sec. 48 212 5.88 
Swing 8 24 , 30 , 10 sec. 48 a ae SL Ne 
Seatac | normal 48 ~-1.64 @.22 
Vibrate 1 20 cps, l 2 2 -4,63 9255 
Sttacies | normal 12 =3 11 ie Zo 
Vibrate 2 20 cps, l g 10 -3.81 Sg) 
Static 1 normal 18 ~3,03 22 
Static 1 normal We -1,69 52,66 
Shock 1 0.3) “fab EZ -2.36 8,45 
Static l normal 60 -2.19 10,06 


Oo O : ~O 
*5°, 10 , 30 sec. represents an aSymmetric swing of 5 off 
vertical in one direction and 10 off vertical in the opposite 
direction with a period of 30 seconds. 
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TABLE IX 
TEST RESULTS FOR 


SULZER MODEL D5 OSCILLATOR 
SERIAL NO, 22 


. 


Test Condition Length Drite Rat¢ Standard 

(hours) Af/F (pp 10 _) Deviation 

hour (x 10) 
Static l normal 65 -6.68 4.97 
Static 2 side 49 -5.14 7.44 
Static 3 upside down 49 -4.31 G2 
Static 4 side 48 -4,26 LORE Z 
Static 5 front 48 -2.90 Od 
Static 6 back 48 -2.64 Teiay 
Static 1 normal 68 -2.26 8.64 
Swing 1 5°, 10°, 30 sec. * 49 2.41 7.85 
Swing 2 a) 10°, 10 sec. 50 -20 5.20 
Swing 3 24_, 30 , 30 sec. 47 -1.76 4.09 
Swing 4 24°, 30°, 10 sec. 47 -1.54 4.08 
Swing 5 55: 10°, 30 sec. 48 -1.42 4.47 
Swing 6 Dare 10 lO sec, 48 “1.32 3 0g 
Swing 7 24°, sO, 30 sec. 48 hcl 3.09 
Swing 8 24°, 30 , 10 sec. 48 -1,08 4.65 
Stacie 1 normal 48 -0.95 5.00 
Vibrate 1 20 cps, 1 g 12 -4,50 7.46 
Static 1 normal IZ -2.74 27u6 
Vibrate 2 ZO°cps.. ig 10 200 2.50 
Seatic | normal 18 -1.99 4.26 
Static 1 normal 17 -0.55 ees 
Shock 1 O23" fall 12 -0.02 7239 
Static l normal 60 -1.14 4.37 


O O Oo 
*5 , 10 , 30 sec, represents an asymmetric swing of 5 off 


vertical in one direction and 10 off vertical in the opposite 
direction with a period of 30 seconds. 
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6. Observations. 

At the onset of the investigation, the investigator had hoped 
to be able to correlate any observed changes in drift rate or frequency 
stability with the conditions of the particular test during which the 
change was noted. A seemingly ideal situation was present with the 
four Sulzer oscillators, three being identical models and the fourth 
an excellent frequency standard, It had been planned to obtain addi- 
tional information from Sulzer Laboratories concerning the orienta- 
tion of the crystal plane within the oscillators. With this informa- 
tion the investigator had hoped to relate the results of the static 
tests with crystal position. 

The manufacturer stated, however, that the position of the crystal 
plane is not fixed, but varies with each individual oscillator. 

Since the only remaining method of determining the crystal orienta- 
tion would have required disturbing the double-oven assemblies, this 
second aspect of the investigation was abandoned. 

The correlation of changes in frequency stability and drift rates 
with test conditions in the case of the Sulzer equipment was hampered 
by the fact that the limited time available for testing required that 
the tests be conducted during the initial warm-up and aging period 
of the oscillators. As the equipment was new, no data record was 
available for use as a behavior reference. These conditions prevented 
the investigator from making the detailed conclusions originally hoped 
for. Instead, general observations of the behavior of the oscillators 


while under test were made, These observations are: 
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Western Electric Frequency Standard 0-76A/U 


Static Tests. Drift rates varied between 0.04 (pp 107) per 
hour and 1.49 (pp 107) per hour, with the lowest value occurring 
while the oscillator was in the normal position. However, each time 
the oscillator was in this position, a different drift rate was 
observed (0.11, 0.10, 0.04). During a portion of the Static Five 
fest with the oscillator tipped forward on its face the drift rate 
was 0.05 (pp ony 

Stability (standard deviation) varied from 0.35 x 10°” to 
ay 72 or with the lowest value occurring while the unit was in a 
normal position. However, as in the case of the drift rates, each 
time the oscillator was returned to its normal position, a different 
value was observed (0.35, 0.32, 0.57). Once again there were instances 
where the stability while the oscillator was in some position other 
than normal was better than that observed for some of the normal 
position test periods, 

The maximum values of drift rate and standard deviation occurred 
during the recovery periods as discussed in Section 5, 

Swing Tests. The drift rate was certainly affected by the 
swinging motion although during the test period from 12 to 14 
September while the oscillator was in the Static One position the 
drift rate was almost as great as the largest observed during 


swinging. The rates varied from 0.133 (pp To to -0.037 


(pp 10°) /hr. 


5) 





Values of standard deviation were generally greater during 
swinging than while at rest, indicating poorer stability under these 
conditions, 

Sulzer Model 2.5 Frequency Standard 

Static Tests. The drift rate gradually decreased as the crystal 
aged and no apparent effects of static position were observed. 

Stability varied between 4,2 x toes and 11.4 x om with the 
greatest value observed during the Static Six test, while the unit 
was tipped on its back. 

Swing Tests. The drift rate continued to decrease gradually 
with further aging of the crystal. In three out of four instances, 
however, the drift rate increased slightly when the swing frequency 
was increased. 


Stability improved as the swing tests progressed, with values 


Wl 11 


of the standard deviation decreasing from 9.3 x 10 to 3.8 x 10 
In three out of four instances, however, the standard deviation in- 
creased when the amplitude of the swing was increased. 

Sulzer Model D5 Oscillators 

Static Tests. As in the case of the Sulzer frequency standard, 
the drift rate for each of the three Model D5 oscillators gradually 
decreased during the aging period. 

Stabilities varied between 5 x Tore and 10 x Ten with no 


common behavior pattern evident. One oscillator exhibited a downward 


trend, one an upward trend, and one remained fairly constant. 
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Swing Tests, The drift rates generally decreased during the 


swinging tests as aging continued. 


del 11 


Stabilities varied from 3 x 10° to15 x10. The lowest 
values of standard deviation for each oscillator were obtained during 
one or another of the swing tests and not during a static period. 
In three out of four instances when the swing frequency was increased, 
the stability of various oscillators improved. The exceptional cases, 
however, did not occur simultaneously. 

The stability of two of the oscillators, serial no. 19 and 
serial no, 22, improved in three out of four instances when the 
amplitude of the swing was increased, 

Vibration and Shock Tests. All drift rates were increased by 


11 
three or four parts per 10 per hour during the initial vibration 


esi) 
11 


e 


Standard deviation varied between 2.2 x 10." and 10.1 x 107 
All standard deviations were increased during the shock test; the 
amounts of increase varied from one to four parts per lOve 

In conclusion, the drift rate and frequency stability of the 
frequency standards were apparently affected to a small degree by 
swinging motion, although this motion had no definite effect on the 
Sulzer Model D5 oscillators. The D5 oscillators were affected 
slightly by shock and vibration. In only a few cases, however, such 
as the change in stability during the shock test, or the change in 


drift rate during the initial vibration test, did all three Model D5 
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oscillators respond in a similar manner. Generally, opposite re- 
actions occurred, i.e., during a given test period, drift rates of 
two oscillators would increase while the third decreased, or the 
stabilities of a different pair of oscillators would decrease while 
the third increased or remained constant. 

In spite of the lack of positive correlation of changes in 
drift rate and stability with conditions of position and motion, 
the results of the investigation were fruitful in that they indicated 


a lack of vulnerability of the oscillator units to motion in general, 
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ip Recommendations for further study. 

Since the initial warm-up and aging period of the oscillators 
is now over, and a backlog of behavior data has been accumulated, 
further investigation should produce more positive results. A repe- 
tition of the tests described in Section 3 would be in order, as 
would be further testing with the Sulzer D5 oscillators mounted in 
shock mounts. 

Test periods should be at least 48 hours in length and, if 
time permits, each test should be preceded and followed by a similar 
amount of time with the oscillators in a normal position, Record- 
ing equipment utilizing 100-division paper should be used to record 
the phase comparator output so that with ten micro-seconds full- 
scale deflection, phase shift at 100 kilocycles may be read to one- 


tenth or one-twentieth of a micro-second. 


Bye, 
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APPENDIX 


A COMPUTER PROGRAM FOR FITTING A STRAIGHT LINE 


TO DATA POINTS AND COMPUTING THE STANDARD DEVIATION 


OF THE POINTS ABOUT THE FITTED LINE 
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104 
999 


998 


100 


996 


30 


40 


997 


PROGRAM LINEFIT 
THIS PROGRAM FINDS THE SLOPE AND Y-INTERCEPT OF THE BEST STRAIGHT 
LINE (IN THE LEAST SQUARES SENSE) FITTED TO A SET OF DATA POINTS. 
THE X DATA ARE ASSUMED TO BE WITHOUT ERROR. TRANSFORMATIONS TO 
EITHER OR BOTH VARIABLES CAN BE MADE. THE PROGRAM WILL ALSO 
COMPUTE THE STANDARD DEVIATION op THE POINTS. 
USAGE 
1+ CALLING SEQUENCE 
MAIN PROGRAM. THE FOLLOWING DATA CARDS ARE READ. 
A. FORMAT (314) - FIRST CARD 
N» RIGHT-JUSTIFIED NUMBER OF: DATA POINTS. 
ICODEX,. TRANSFORMATION CODES FOR X 
ICODEY, TRANSFORMATION CODES FOR Y 
B. FORMAT (12F6.0) - SUBSEQUENT CARDS. 
X(1),Y¥(1)2X€2),Y(2),5 ETC. 
C. TRANSFORMATION CODES 
=1 NONE : 
=2 VARIABLE REPLACED By SQUARE: ROOT OF VARIABLE 
©3 VARIABLE REPLACED BY SQUARE: OF VARIABLE: 
=4 VARIABLE REPLACED BY LOG TO BASE 10 OF VARIABLE 
=5 VARIABLE REPLACED By LOG:-TO BASE E OF VARIABLE 
41111 ; 


TEST PERIOD FROM 1700 4 OCT TO 1500 7 OCT, S1, STATIC ONE, FINAL 


DIMENSION X(200),XSAVE(200),Y(200),YSAVE(200),YFIT(200), DIFFS(200 
READ 999, Ns ICODEX, ICODEY 

FORMAT(314) 

PRINT 998,N> ICODEX» ICODEY 

FORMAT(4H1N =14,3X8HICODEX =I2, 3x 8HICODEY =I12//) 

READ 100. (X(1),YC(1),124,N) 7 Se 
FORMAT(12F6.0) “Sts 

DO 996 I=4, 
XSAVE(])= : 
YSAVE(I)= 
YSUM 
SUMX 
SUMY 
S U M x Y = 0 @ 

SUMXX = 0. 

CALL TRANS(N,X,Y,1CODEX, ICODEY) 

XN = N 

DO $0 | = 1+ +N 

SUMX=SUMX*#X( 1) 

SUMY = SUMY+Y(1[) | 
SUMXX = SUMXX + (XCD) #X¢C1)) = 
SUMXY=SUMXY+ (XCL)#*YCT)) 

DENOM = XN*SUMXX=SUMX#SUMX 

ANUM = SUMY * SUMXX - SUMX * SUMXY 

A = ANUM / DENOM | 
BNUM= XN®SUMXY=SUMX*SUMY eves: 

8 = BNUM / DENOM a. 

DO 40 J. = 1.N 

YC = A+B * X(J) 


N . . ' ‘ 
1) “yy 
(1) 


oo ou 
Lon | 


-DIFF = Y(U) - YC 


YSUM = YSUM + DIFF eDIFF 


S22 = YSUM / (XN - 1.) 


S1G=S2/DENOM 

SRA=SQRTF (SUMXX*#SI1G) 

SRB=SQRTF (CXN*SIG) : 

PRINT 997, As B, SRA,- SRB 

FORMAT(24HQ STRAIGHT LINE’ FITTING// 25H USING THE FORM Y=BX + A, 


1/ 4H A = E14.8/ 4H B = E14.8/ 26H STANDARD DEVIATION OF A 8: 644.8 
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90 
995 
99 


1000 


106 


401 
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2/ 26H STANDARD DEVIATION OF B = E14.8) 


SUMFIT=.0 
SD WPRrS=0. 
SDIFFSQ=0. 


DO 90 [=1,N 

YFITCI)=BeXC 1) +A 

DIFFSCI)SYFITCI)-YC1) 

SDIFFS=SDIFFS+DIFFS(]) 

SDIFFSQ=SDIFFSQ+DIFFS(I) «#2 

SUMFIT=SUMFIT+ DIFFS(1) 

PRINT 99,SUMFIT,(X( I), YSAVECID,YCID-YFITCI) DIPFSCI), ledtsNn> 
FORMAT(QHOSUMF IT =E12¢65///10X1THX1I4X1HY9OXLOHYTRANSPOSESBX4HYP IT11X5) 


LDIFFS///(5E15.9)). 


SIGMA=SQRTFCCSDIFFSQ-SDIFFS#*2/XN)/(XN -1.)) 
PRINT 1000-SIGMA 


FORMAT (1HO»7HSIGMA =E15.8) 

GO TO 194 

7 OP pa 
END 


SUBROUTINE TRANS(N, X, Ys ICODEX, ICODEY) 
DIMENSION X(200), Y(200) 

GO' TO ( 105 20+ SO» 40+ 59), ICODEX 

DO 21 I = 15N 


X(T) =SORTFCABSFC(X(1))) 


GO TO 10 

DO $1 1 = 1.N 

XC(I) = XCT) & X¢J) 

GO: TO 10 

DO 41 | = 1,N 

X(T YSLOGLOF(CX(1)) 

GO TO 10 

DO 51 I = 15N > 

X( I) =LOGF(XC(I)) 

GO TO (100, 200, 300, 400, 500), ICODEY 
DO 201 | = 1,N 
YCI)=SQRTFCABSFC(YC(I))) 
RETURN 

DO 5o1 I = 1,N 

YC(I) = YC) ® YCT) 
RETURN 

DO 401 I = 1,N 
Y(ITVSLOGLOFCYC(T)) 
RETURN 

DO 501 I = 1,N 


Y(1)=LOGF(Y(T)) 
END 
END . 
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N= 68 ICODEX = 1 


STRAIGHT LINE FITTING. 


USING THE FORM Y=BX + A, 
A =-.99372827E+93 
B =-.239753526+91 
STANDARD DEVIATION OF A 
STANDARD DEVIATION OF 8B 


te 


ICODEY 


SUMFIT = .40431E-06 
x Y 
-10000E*01 -.90400E+035 
-20000E*01 -.91790E+05 
~S0000E*01 -.91000E+03 
-40000E*01 -.91000E+03 
-20000E*01 -.92490E+938 
-6O000E+01 -.91700E+03 
~-70000E+01 -.91000E+038 
~-80000E+01 -.93100E+03 
-ID000E+01 -.95100E+03 
-10000E*02 -.93100E+03 
-11000E+02 -.93100E+03 
-12000E+02 -.93100E+03 
-13900E+02 -.93100E+03 
-14000E+02 -.935100E+03 
-15000E+02 - 938Q0E+03 
-16000E+02 -.944990E+03 
-17000E+02 -.94490E+03 
-18Q900E+02 -.94490E+03 
-19900E+02 -.94490E+03 
~20000E+02 -.94490E+05 
-21000E+902 = .95890E+93 
»-22000E+02 -.96500E+03 
-23000E+02 -.96500E+03 
-294000E+92 -,. 95890E+03 
~25000E+02 -. 96599E+93 
-26Q000E+02 -,.95890E+03 
-27000E+02 -,.96590E+03 
~-28900E+02 -,.96590E+038 
~-29000E+02. -.97200E+038 
~S0000E+02 -.96590E+03 
~S1000E+02 -.97900E+03 
~S2000E+02 -.97200E+03 
~S3000E+02 -~.97900E+038 
~54000E+02 -.96500E+03 
-355000E+02 -.97900E+03 
.5S6000E+02 -.99300E+03 
~-37000E+02 -. 98600E+03 
~-SB8000E+02 -.99300E+03 
-59000E+02 -.98600E+03 
~-40000E+02 -.99300E+03 
~-41000E+92 - 993Q0E+903 
~-42900E+02 -.993SQQ9E+Q93 


Cad 


YTRANSPOSE 


= .904Q90E+Q5: 
7 .91700E+93: 
= .91000E+035. 
‘=,91000E+03 
=, 9249 Q9E+05: 
-.91700E+93: 


-,.91000E+05 


‘* .95100E+95: 
-=,95100E+03: 
-.935100E+03. 
-.9S100E+03. 
-,93S100E+93. 
~,93100E+ 93: 
~.95100E+93: 
~ 9S8Q0E+QS5: 


-.944900E+903 


-,94490E+ 93. 
-.94499E+Q93. 
= ,944Q9Q9E+93. 
-,94499E+93:- 
‘= 958QQ9E+93 
= ,96500E+03: 
= ,965Q99E+93: 


-.95890E+05 


=~ ,96599E+93. 
--, 958Q0E+93 


-,965Q90E+93 


-=,965Q99E+93. 
-,.97200E+ 93. 
-, J65Q90E+ 93: 
- .97900E+903~- 


~.97200E+038 


-.97900E+03. 
= ,96500E+93. 
-.979Q00E+903. 
-,.99300E+Q93. 
-. 98600E+03 
-.99500E+03 


-. 9B600E+93 


‘= 99SQQ0FE+035: 


- .993SQ0E+93 


‘™, 99SQ0E+Q3. 
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-,.99694E+93 


"= .99834E+93 
- *.91965E+93. 


~.91296E+93 
-.91527E+93 
~.91757E+93 
-.91988E+93 
~ 92219E+Q93 


-.92459E+03 
~.92680E+03. - 


-.92911E+903 
-.95142E+93 


- ISS7SE+93 


as 9S60SE+98 
-,. 938354E+93 


-.94965E+903. 


-.94296E+93 


-.94526E+93. 


-.949757E493 
-,9498B8E+Q93 
- .9S219E +93 


a -95449E+93. 


-.95689E+03 
-.95911E+93 


~,96142E+ 93. 


-,.96372E+93 
-.966903E+93 


-.96854E+93. 
7 e97065E+08 


' ' = ,97526E+03 


-.97757E+03 
-.9/7988E+ 93 
-.98218E+03 
-,.98449E+93 


-.98680E+05. 
aeags ~9B9OLLE+03. 
= ©991491E+903. 


-,99372E+93 
-,.9960SE+93 
-~.99BS4E+93 


“= .100Q05E+04 


DIFFS 


—2QOS5BE+91 
~86567E+91 
-65988E+99 
—29584E+91 
8735416401 
~D754BE+90 
»98810E+91 
»881149E+91 
~650359E+91 
»-41964E+91 
~18888E+91 
~41869E+90 
~27262E+91 
~00S38E+91 
»54130E+00 
05S5512E+91 
»10456E+01 
-12639E+91 
as ~55/714E+91 
28/799E+91 
58135E+91 
°10506E+92 
~81984E+91 
-~.11091E+91 
~SSBSSE+91 
-,97242E+91 
-.10517E+91 
~ SSSOSE+91 

»15532E6+91 
~,79543E+91 

~S7SB1LE+91 
-,55694E+91 
-.87693E+90 
-.17184E+92 
- 54920E+01 
»62005E+01 
-.51071E+901 

°15854E+91 
~./7221E+91 
- , S0297E+91 
~,ISS72E+914 
-,76447E+04 


5 





~43000E+02 -.10000E+04 
»44900E+02 -.10000E+04 
~45000E+02. -.10070E+04 
-46000E+02 -.10140E+04 
»-47000E+02 -.10070E+94 
»48000E+02 -.10140E+04 
~49000E+02 -.10280E+04 
-20000E+02 -.10280E+04 
-21000E+02 -.10210E+04 
52000E+02 -.10280E+04 
-°3000E+02 -.10210E+04 
~-24000E+02 -.10210E+04 
~-25000E+02 -.10350E+04 
56000E+02 -.1903506+04 
7 D00EtO2 -.10420E+04 
-58000E+02 -.10420E+04 
59000E+02 -.10420E+04 
-60000E+*02 -.10420E+04 
-61000E*02 -.10420E+04 
»-62000E*02 -.10420E+04 
.63000E+02 -.104920E+04 
-64900E+02 -.10490E+04 
-65900E+02 -.10560E+04 
~66990E+92 -.19930E+04 
.67900E+*02 -.105690E+04 
-68900E+02 -.10659E+04 
SIGMA = 


-56589160E+01 


-. i0000E+04 


-.10000E+04: 


-.10070E+04 


‘--,10140E+04 


-.10070E+04 
~,10140E+04 


—-,10280E+04 
-=,10280E+04 
‘=,10210E+04 

-,10280E+04: 


~,10210E+04 


-,10210E+04: 


-.10550E+04 


=,10550E+04 
-.10420E+04: 


-,.10420E+04 
-~.10429E+94 


=.,10420E+94 
=.104290E+ 94 
--,190429E+94 
‘= ,10420E+04 
‘7 .10499E+9 4: 
= .10960E+04: 
-.19939E+94: 
-=,10690E+04: 
°,10639E+04: 
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-,100350E+04 
- .10093E+ 94 
-.10076E+04 
-.10099E+04 
-.10122E+94 
-.10145E+04 
-.10168E+04 
-.10191E+04 
-.10214E+04 
~.10257E+04 


1.10260E+04 © 


-.10283E+04 
-.10306E+04 
-.10330E+04 
-~.10553E+04 
~.10576E+04 


-.10599E+94: 
--"=.10422E+94 
-°.10445E+04 
= .19468E+94 


-,10491E+04 


: ~.10914E+94 
‘ -.109357E+094 
#.190569E+94: 


-,1090585E+94 


-.10696E+94: 


-.29925E+01 
- ,52598E+91 
»96756E+00 
~41251E+91 
,51824E+91 
~48996E+00 
.11203E+92 
»8B950E+01 
~41257FE+900 
42799E+91 
~D20276E+01 
~7SS52E+91 
043575E4+91 
»20498E+01 
67422E+01 
—, 44347E+01 

6214272E+91 
~,18938E+00 
29879E+91 
eF#7955E4+91 
eALOSOE+01 
291956491 
»22819E +91 
»69744E491 
610667E+92 
»2859SE+Q1 























